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The dynamics of a doped hole in the orthogonal-dimer spin system is investigated system- 
atically in one, two and three dimensions. By combining the bond-operator method with the 
self-consistent Born approximation, we argue that a dispersive quasi-particle state in the dimer 
phase is well defined even for quasi-two-dimensional systems. On the other hand, a doped hole 
in the plaquette-singlet phase hardly itinerates, forming an almost localized mode. We further 
clarify that although the quasi-particle weight in the dimer phase is decreased in the presence 
of the interchain coupling, it is not suppressed but even enhanced upon the introduction of the 
interlayer coupling. 
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Recently, low-dimensional quantum spin systems with 
frustration have attracted much interest. Among 
them, a quasi-two-dimensional (2D) spin gap compound 
SrCu2(B03)2 1:2 ' is outstanding in its unique features. 
The crystal structure of the compound is such that Cu 2+ 
ions with spin S = h are located on the orthogonal- 
dimer lattice, resulting in the characteristic magnetic 
properties. 3 ~ 7 ) Miyahara and Ueda 8 ) pointed out that 
the orthogonal-dimer lattice is equivalent to the square 
lattice with some diagonal bonds, which is referred to as 
the Shastry-Sutherland model. 9 ) (See Fig. 1) The sys- 
tem has been providing a variety of interesting issues, 10 ) 
e.g. quantum phase transitions, 8, n ~ 19 ) unique triplet ex- 
citations, 12, 20 ~ 25 ) magnetization plateaus. 14, 20 ~ 22,26 ~ 29 ) 
It was found that the ground state of the compound 
SrCu2 (1303)2 is in the dimer phase close to the quan- 
tum phase transition point. 8 ) Therefore, the application 
of the magnetic field 1,3 -* or the pressure, 7 - 1 could trigger 
a quantum phase transition. Also, some reports claimed 
that the compound SrCii2 (1303)2 has a rather large inter- 
layer interaction to explain the experimental data such 
as the susceptibility and the specific heat. 12,30 -* 

The study of a single hole doped into the above 
orthogonal-dimer system is important, 31 ) since its dy- 
namics determines the profile of the photoemission spec- 
trum. In comparison with the extensive theoretical inves- 
tigations of static quantities, such one-particle spectral 
properties have not been studied systematically, except 
for Vojta's work on the dimer phase of the 2D Shastry- 
Sutherland model. 31 ) In particular, it is desirable to 
clarify how the interlayer coupling, which turns out to 
be important for the static quantities, affects the pho- 
toemission spectrum. Moreover, such investigations of 
a doped hole may provide a starting point to study the 
system with finite density of holes, which may be realized 
by the chemical substitution of some elements. Although 
simple substitution of Zn ions for Cu ions results in lo- 
calized magnetic moments, 32 ) it is an interesting issue to 
realize a metallic orthogonal-dimer system experimen- 
tally. 



Motivated by the above topics, we investigate the 
dynamics of a doped hole in the orthogonal-dimer 
system systematically in one, two and three dimen- 
sions. By means of the self-consistent Born approxi- 
mation (SCBA) 33 ~ 38 ) with the bond-operator represen- 
tation 31, 37 ~ 41 ) and also the exact diagonalization (ED), 
we discuss how the characteristic orthogonal-dimer struc- 
ture affects the motion of a doped hole. 

This paper is organized as follows. In Sec. 2, we in- 
troduce the orthogonal-dimer spin system, which is rel- 
evant to describe the system of SrCu2 (603)2- In Sec. 
3, we briefly summarize the bond operator representa- 
tion for quantum spins, which is combined with the self- 
consistent Born approximation to treat a doped hole. 
This method is used for the analysis of the dimer phase, 
whereas the ED method is also employed to treat the hole 
dynamics in the plaquette phase. In Sec. 4, we study the 
dynamics of a doped hole systematically in one, two and 
three dimensional orthogonal-dimer systems, and clarify 
the role of the lattice geometry on the dynamics of a 
doped hole. A brief summary is given in Sec. 5. 

2. model 

We investigate the orthogonal-dimer spin model, 12, 42, 43 ) 
which is schematically drawn in Fig. 1. To deal with the 
system with a hole, we consider the following t-J Hamil- 
tonian, 
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Tk,-o)]^\ where cSt2 annihilates 



(creates) an electron and nj j(T is the electron number 
with spin ct(=T,4) at the ith site. The corresponding 
spin operator is defined by S, = \ J2 a ,a> A,o T a,a lC i,a' ■> 
where r is the Pauli matrix. We introduce four types 
of hopping integrals t K > and exchange couplings 
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Fig. 1. Orthogonal-dimer system for SrCu2(B03)2, which has al- 
ternating orthogonal-dimers in all directions along the a, b and 
c axes: (a) the 2D Shastry-Suthcrland model in the a-b plane. 
Square unit cell containing four spins is shown as a long dashed 
line. If we set J3 = tz = (dotted line), the system is reduced 
to an orthogonal-dimer chain with alternating dimcr-plaqucttc 
structure; (b) configuration of the dimers along the c-axis, which 
forms another type of orthogonal-dimer chain. This is equiva- 
lent to a ladder model with additional exchange interactions for 
diagonal bonds, which is referred to as tetrahedra chain. 



J K > 0(k — 1,2,3,4), which are specified by the bond 
pairs D K = {J K ,t K } shown in Fig. 1. Note that this 
generalized model has the orthogonal-dimer structure in 
all three directions, thereby allowing us to discuss the 
hole doping effect systematically in one, two and three 
dimensions. When J2 = J3 and J4 = 0, the system is re- 
duced to the well-known 2D Shastry-Sutherland model, 9 > 
static properties of which have been intensively studied 
so far. 8 ' n ~ 24 ' 26-28) rpj^g gyg^gj^ j g f ur ther reduced to a ID 
orthogonal-dimer system if we consider the chain system 
along the a or b axis (lower panel of Fig. 1(a)). Also, the 
ID system along the c axis forms another orthogonal- 
dimer chain, as shown in Fig. 1(b). Magnetic prop- 
erties of these chain models have been studied in de- 
tail 44 ~ 47 ' 49 ' 50 ' 5 °- 55 ) 

In contrast to well-studied static properties, system- 
atic investigation of one-particle spectral properties 31 ) is 
still lacking, which we will address in the following sec- 
tions. 

3. Bond operator approach 

In this section, we introduce a method to treat the dy- 
namics of a doped hole. Since we are mainly concerned 
with the dimer phase, it is convenient to use the bond 
operator representation, 31, 37 ~ 41 ) which may describe the 
ground state as well as low-energy excitations in the 
dimer phase rather well. The essence of this method 
is to represent the Hamiltonian in terms of " bond opera- 
tors" , which act on each dimer pair formed by J\ , instead 



of original operators for spins and electrons. Since the 
generalized Shastry-Sutherland model has four sites in a 
unit cell, we introduce the bond operators defined as, 
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for a horizontal bond at the n-th cluster in Fig. 1. 
The operators s(s'),t a (p a )(a — x, y, z) obey the bosonic 
commutation relation, while the operator a the fermionic 
anticommutation relation. The notations u,d,l,r spec- 
ify the up, down, left and right sites in the unit cell as 
shown in the lower panel of Fig. 1(a). The ket states 
are labeled by the configuration of electrons sitting at 
up (left) and down (right) sites of the vertical (horizon- 
tal) bond [o represents an unoccupied (hole) site]. To 
restrict the Hilbert space to the physically sensible one, 
we introduce the following constraints, 

4 s " + J2 *i,„*a,n + J2 a i,n,a a i,n,<? = 1. (6) 
ct a,i—u.d 

s'U'n + Y,PknP<*,n + 4,n,cr a i,n,<T = 1- (7) 

We are now considering the dimer phase for the host 
spin system, which is realized by the condensation of the 
singlet-bond operators, s n and s' n . The remaining triplet 
operators describe spin excitations in the host system, 
giving Hj in the Fourier space, 

Hj = ~J X N 

+^2{^ikt f ak t aM +uj2kpi M p a ,k} (8) 
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with 



parameters is given as, 



*t . = </Xw 

k 



(9) 



where uj\k — ^2k = Jl, N is the number of unit cells and 
k is the wave vector. When we have derived the Hamil- 
tonian (8), the interactions among triplets have been dis- 
carded due to their minor importance. We note that the 
triplet excitations uj\ and u>2 obtained do not depend on 
the wave number, implying that the orthogonal-dimer 
structure forbids the motion of the triplet. In fact, it is 
known that the triplet excitation does not have the dis- 
persion up to the fifth order in the inter-dimer coupling 
for the 2D model. 8 -' To make our discussions more quan- 
titative, we use the triplet excitation energy obtained by 
the perturbation expansion up to the fourth order, 
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We now consider the hole dynamics in the orthogonal- 
dimer system. To this end, we define the retarded Green 
function for physical fermions in the 4x4 matrix form, 

G a (k,i) = -iQ(t)(D\{^ t<T (t)^l a }\D), (12) 

where 4 )(T = (4,k )CT . 4,k, CT . 4,k, CT > 4,k )CT ) and \ D ) is the 
dimer singlet ground state. 8 ' We employ here the SCBA 
to obtain the self-energy of the Green functions. It 
is known that this approximation works rather well as 
far as a single hole doped in quantum spin systems is 
concerned. 33 ~ 38 ) Following standard procedures in the 
SCBA, we compute the spectral function of the physical 
fcrmion defined by 



A(k,u) = Im Tr G(k,u), 

7T 



(13) 



which directly gives the photoemission spectrum. This 
function satisfies the sum rule, J A(k, oj)duj = 2. 

This completes the formulation of the bond-operator 
treatment of the dimer phase. We will also employ the 
ED calculation of the finite-size system to check the va- 
lidity of the bond-operator approach. Moreover, we will 
see that the latter ED calculation is efficient to discuss 
the hole dynamics in another singlet phase of the model, 
i.e. the plaquette phase. 

Here we make a comment on the exchange couplings 
of our generalized model. If we consider the t-J model 
(1) is to be derived from a strong-coupling limit of the 
Hubbard model with on-site repulsion U , the ratio of the 
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We use these relations when the ratio of the exchange 
coupling and hopping is altered. 

4. Hole dynamics in orthogonal-dimer system 

4-1 orthogonal-dimer chain 

We begin with the ID spin chain ( J3 = J4 = 0) shown 
in the lower panel of Fig. 1(a), which may have es- 
sential properties inherent in the orthogonal-dimer sys- 
tems. 44 - 47 ' 50 ) Before presenting the computed results, 
we note that the ground state of the chain model be- 
longs to the dimer (plaquette) singlet phase for J2/J1 < 
J c {> Jc) with J c = 0.819. 44 - 47 ) Shown in Fig. 2 is the 
spectral function A(k, u>) for the dimer phase, which has 
been calculated by combining the bond operator method 
with the SCBA. Since there are four sites in the unit cell, 
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Fig. 2. Spectral function A(k, ui) for the orthogonal-dimer chain 
with Ji = 5, J2 = 5 X 0.635 and J3 = J4 = (dimer phase). 
The contribution of the dispersionless mode is shown by the 
dashed line. 

we have four kinds of clear peaks in the spectral function: 
three dispersive modes and one dispersionless mode. We 
refer to the lowest mode having clear dispersion as the 
quasi-particle state. We thus say that the quasi-particle 
state is well defined and is mobile rather freely in the 
background of the dimer singlets. 

In the plaquette phase, the spectral function A(k,u>) 
shows a similar four-peak structure, as seen in Fig. 3, 
which is obtained by the ED calculation of the finite sys- 
tem (Ji = 5, J2 = 5 x 0.85). In this case, however, the 
lowest quasi-particle peak is dispersionless, in contrast 
to the dimer phase, implying that a doped hole is for- 
bidden to itinerate in the chain system. This results in 
a completely localized quasi-particle state. 

We now clarify why the above characteristic features 
appear in the spectral function, by taking into account 
symmetry property inherent in the orthogonal-dimer 
structure. This consideration may be also helpful to dis- 
cuss the 2D and 3D cases. Notice first that the sys- 
tem has symmetry of space inversion P for each vertical 
dimer-pair; namely if the pair forms a singlet (triplet), 
we have the eigenvalue P = 1 (P = — 1). Then, the 
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Fig. 3. Spectral function A(k, ui) for the chain in the plaquette 
phase (Ji = 5, J2 = 5 X 0.85) obtained by the ED for the small 
system (16 sites). 



Hilbert space of the Hamiltonian is classified into each 
subspace specified by a set of the eigenvalues [{Pi}]- In 
the dimer (plaquette) phase for the undoped case, the 
singlet (triplet) state is realized on every vertical bond. 47 ) 
Therefore the ground state of these phases is respectively 
given by the lowest state in the subspaces [111 • • •] and 
[III • • • ], as shown in Fig. 4(a). 

Now consider the hole doping effect. In the dimer 
phase, a certain dimer singlet may be broken by a doped 
hole, but the bonding state (symmetric superposition) 
is formed between the doped hole and the correspond- 
ing unpaired spin. Therefore the ground state in the 
doped case is still in the same subspace [111 • • •], since 
the bonding state on a vertical bond has the eigenvalue 
P = 1 ( Fig. 4(b)). This symmetry property renders the 
hole state mobile, making a well defined quasi-particle 
state. Two other dispersive modes in Fig. 2 belong to 
this class of P = 1 . On the other hand, a dispersionless 
hole state is realized by creating an anti-bonding state 
(P = —1) at a broken-dimer site, which thus belongs 
to the subspace [1111 • • •]. Since the different inversion 
symmetry between neighboring pairs prohibits the mo- 
tion of a hole, leading to a completely-localized excited 
state. 

Similar consideration can be applied to the hole dop- 
ing in the plaquette phase. In this case, however, the 
energy level of the anti-bonding hole state (P = — 1) is 
higher than that for the bonding state (P = 1). There- 
fore, the ground state is not realized in the subspace of 
[1111 • • • ] but in [1111 • • • ], as indeed seen in Fig. 4(b). 
Therefore, we end up with a completely-localized quasi- 
particle (lowest-energy) state upon hole doping in the 
plaquette phase. 

4-2 2D Shastry- Sutherland model 

We now turn to the 2D case. In the following, we dis- 
cuss how the spectral function of a doped hole is changed 
when the interchain coupling J 3 is introduced. We be- 
gin with a quasi-particle state in the dimer phase. The 
spectral function obtained by the SCBA is shown in Fig. 
5. It is seen that the introduction of the interchain cou- 
pling increases the band width of the lowest mode while 
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Fig. 4. The energy of the lowest state in the distinct subspaces (a) 
undoped case, (b) one-hole case. The data are obtained by the 
ED for the system (24 sites) with periodic boundary conditions. 



keeping its dispersive peak. This implies that the quasi- 
particle state is well defined even in the 2D Shastry- 
Sutherland model, in accordance with the results of Vo- 
jta. 31 ) On the other hand, other peaks in the high energy 
region are considerably smeared with increase of the in- 
terchain coupling. Such smearing effects are induced by 
scattering by triplet bosons. 

Now let us turn to the plaquette phase. Note that 
this type of plaquette-singlet phase may exist even in the 
2D Shastry-Sutherland model due to strong frustration, 
which is believed to be sandwiched by the dimer phase 
and the antiferromagnetically ordered phase. 11, 18, 19 ) We 
discuss here the dynamical properties of the doped hole 
in the plaquette phase. By using the ED calculation of 
the 4x4 small cluster, we obtain the spectral function 
as shown in Fig. 6. The introduction of the interchain 
coupling J3 has little effect on the quasi-particle state 
(lowest-energy state) except for the simple energy shift. 
In particular, it should be noticed that the dispersionless 
quasi-particle state can persist even in the 2D system. 
We also note that higher dispersive peaks are not so much 
smeared in comparison with those in the dimer phase. 
This is partially due to the finite-size calculation we have 
used for the plaquette phase. 

For reference, we plot the band width of the quasi- 
particle state in Fig. 7. In the dimer phase, the band 
width is roughly proportional to hopping t, and thus it 
increases in proportion to the square root of J2/J1 since 
the ratio of hopping and exchange coupling is fixed as 
in eq. (14). At the critical point J c , the band width 
shows discontinuity, reflecting the first-order quantum 
phase transition. Beyond the critical point J c , the doped 
hole hardly itinerates, as mentioned above. In particular, 
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Fig. 5. Spectral function A(k, ui) in the dimer phase (Ji = 5 and 
J 2 = J x x 0.635). 
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Fig. 6. Spectral function A{k,u>) in the plaquette phase. The 
choice of the parameters are Ji = 5 and J2 = J3 = Ji X 0.85, for 
which the plaquette singlet phase should be stabilized. 11, 18 > 19 ) 



in the chain case (J3 = 0), it is completely localized. 
By increasing the interchain coupling J3 , the doped hole 
can itinerate in the system even for the plaquette phase, 
giving rise to a slight increase of the band width. 

Summarizing the above results, we can say that the 
one-particle (photo-emission) spectrum provides an ef- 
ficient probe to distinguish the two spin-gap phases 
clearly; i.e. a dispersive (almost localized) quasi-particle 
state characterizes the dimer (plaquette) phase. We wish 
to note that similar behavior appears in a triplet spin ex- 
citation spectrum in the undoped case, for which char- 
acteristic properties in the two phases are interchanged. 
Namely, a triplet excitation forms a completely local- 
ized (dispersive) mode for the dimer (plaquette) phase 
in the ID case. Even in 2D, it has been shown that a 
triplet excitation in the dimer phase is almost localized. 8 ) 
Therefore, we can efficiently use these properties to an- 
alyze the photoemission spectrum and/or the neutron 
scattering experiments to figure out whether a plaquette 
phase can be really realized experimentally by changing 
the pressure, the magnetic field, etc. 

4-3 Effects of the interlayer coupling 

Finally, we discuss to what extent the above char- 
acteristic properties can persist in the presence of the 
inter-layer coupling. Some recent reports claimed that 
the compound SrCu 2 (B0 3 ) 2 has rather large interlayer 
couplings, which may not be negligible to explain the ex- 
perimental findings. 12 > 30 ) Nevertheless, it is known that 
some of physical quantities can be explained rather well 
by the 2D model. We address this question on the hole 
dynamics in the quasi-2D orthogonal-dimer system. We 
are mainly concerned with the dimer phase below in this 



section. 

We show the computed spectral function of a doped 
hole in the quasi-2D case in Fig. 8, where we have chosen 
two different interlayer couplings J4/J1 = 0.09 and 0.2 
with fixed J2/J1 = J3/J1 = 0.635. A remarkable point is 
that the introduction of the interlayer couplings J 4 and t\ 
does not obscure the spectral function, and seems to even 
sharpen its characteristic peak structure in some energy 
region. We show below that this remarkable behavior 
is caused by the specific arrangement of dimers in the 
direction along the c-axis. 

To make this point clear, we exploit a simplified model 
of the two-leg ladder with diagonal bonds, which is the 
basic structure along the c-direction (see Fig. 1 (b)). 
Magnetic properties of this ladder model have been al- 
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Fig. 7. Band-width of the lowest energy band (quasi-particle 
state) as a function of J2/J1 with J3/J2 being fixed. 
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Fig. 8. Spectral function A(k, u) in the dimer phase (Ji = 5 and 
J2 = J3 = Ji X 0.635) with the interlayer coupling J4 = J\ xO.09 
and J4 = Ji X 0.2. Both choices belong to the dimer phase. 



ready studied in detail. 49 ~ 55 ' To analyze the hole dy- 
namics in this ladder system, it is convenient to introduce 
the fermionic operators for bonding and antibonding hole 
states, which are defined as 



(15) 



Here the fermionic operators a 1(T ,n and a%,a,n denote hole 
states created at two distinct sites of the n-th rung in 
the c-direction. These two states span two independent 
Hilbert spaces when the hopping and the exchange cou- 
plings are introduced along the c-direction. The Hamil- 
tonian for the bonding hole state is immediately diago- 
nalized, and the corresponding eigenstate and the eigen- 
value are given as, 



|<\,„ T ./,-> = -y=^2exp(ikn)al i(T>n s n \D), (16) 

(17) 



'N' 

it, 

t\ — 2^4 cos k, 



where N' is the number of unit cells. An important point 
is that the doped hole in a bonding (symmetric) state 
can move in the ladder without producing any triplet ex- 
citations. This is a strong constraint coming from lo- 
cal symmetry due to the unique dimer alignment in the 
c-direction. Therefore, the quasi-particle peak appears 
clearly in Fig. 9. On the other hand, an anti-bonding 
(anti-symmetric) hole state a a is completely localized, 
as shown in Fig. 9. We have also confirmed this conclu- 



sion by calculating the spectral function by the ED for 
the small cluster N' = 8, which is shown in Fig. 9 by 
the dashed line. It is seen that the results show good 
agreement with those obtained by the SCBA. 




Fig. 9. Spectral function of tetrahedra chain. (a)Ji = 2 and J4 = 
Jl X0.09 (b)Ji = 2 and J4 = Ji xO.2. The dashed line represents 
the results obtained by the ED. 

The above analysis is valid for the ladder system, but 
most of essential features may appear even in the quasi- 
2D case along the c-direction although the profile of the 
spectrum may be somewhat obscured. This naturally 
explains why the quasi-particle peak in the low energy 
region in the orthogonal-dimer system is not smeared by 
the interlayer coupling. 

To conclude this section, we check how well the quasi- 
particle state is stabilized in 3D by estimating the quasi- 
particle weight defined by 



Z(k) 



1 - Re 



dS(k,u;) 
doj 



(18) 



where £k is the energy of the quasi-particle. Shown in 
Fig. 10 is the weight computed at the T point when the 
ratio of the exchange coupling and the hopping is varied. 
In the parameter region shown in the figure, there are 
two characteristic features, (i) the change in character 
from ID to 2D: By comparing the data of J3 = J4 = 
(dotted line) with those of J4 — (broken line) , wc notice 
that the increase of the interchain coupling results in 
the decrease of the quasi-particle weight in accordance 
with the discussions given above, (ii) Change from 2D 
to 3D: On the other hand, the increase of the inter-layer 
coupling does not lead to the decrease of the weight, but 
has a tendency to increase it, as seen from the solid and 
dash-dotted lines in Fig. 10. This property is related to 
the fact that the motion of a hole in the c-direction is 
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Fig. 10. Renormalization factor at the T point with (a)J2/Ji 
J 3 /Ji = 0.2 and (b)J 2 /Ji = J3/J1 = 0.635. 



quite free, rendering the quasi-particle state more stable. 

5. Summary 

We have investigated the hole dynamics in the 
orthogonal-dimer model proposed for the quasi-2D com- 
pound SrCu2(BC>3)2. In particular, the characteristic 
properties present in the one-particle spectrum have been 
clarified in detail, by systematically studying the one, 
two and three dimensional systems. We have demon- 
strated that the quasi-particle state of a doped hole is 
well stabilized due to the specific orthogonal-dimer struc- 
ture. Furthermore, it has been found that the one- 
particle spectral function provides an efficient way to 
distinguish the two spin-gap phases; i.e. a dispersive 
(almost localized) quasi-particle mode characterizes the 
dimer (plaquette) phase. This remarkable property is in- 
herent in the orthogonal-dimer system, which stems from 
a strong constraint due to local symmetry. 

We expect that the plaquette singlet phase could be 
realized experimentally by applying the pressure or the 
chemical pressure with substitution of certain elements. 
In order to judge which type of spin-gap phase is realized 
in such case, an efficient probe is necessary to distinguish 
them properly. We have found here that the photoemis- 
sion spectrum can be an efficient method to specify the 
nature of the spin-gap phases. 
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